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Abstract. We show that the fluctuations of the partial current in two dimensional diffusive sys- 
tems are dominated by vortices leading to a different scaling from the one predicted by the hy- 
drodynamic large deviation theory. This is supported by exact computations of the variance of 
partial current fluctuations for the symmetric simple exclusion process on general graphs. On a 
two-dimensional torus, our exact expressions are compared to the results of numerical simulations. 
They confirm the logarithmic dependence on the system size of the fluctuations of the partial flux. 
The impact of the vortices on the validity of the fiuctuation relation for partial currents is also 
discussed in an Appendix. 

1. Introduction 

Recently, it has been shown how to compute the large deviation function of the current in 
one dimensional diffusive systems [2]- [2]. The hydrodynamic large deviation theory [21 \T9\ IT2] . 
yields explicit expressions for the large deviation function as well as the cumulants of the current 
fluctuations (under some stability condition [5l [7]). The same hydrodynamical approach applies 
in principle also to currents in higher dimension. In the present paper we show however that this 
approach does not always catch the correct scaling of the large deviations or of the cumulants of the 
current in higher dimensions. This will be made explicit in the case of the 2 dimensional symmetric 
simple exclusion process (SSEP). 

For a one dimensional diffusive system of length L in contact at its left end with a reservoir at 
density pa and at its right end with a reservoir at density pb, one can consider the total net number 
Q(r) of particles leaving the left reservoir during a time interval r. This number Q{t) fluctuates 
in time and one expects that in the long time limit 



Pro 



~ exp [ - rGL(J;Pa,P6)] (1-1) 



where Gl{J] PaiPb) is the large deviation function of the flux through the system. In fact Gl does 
not depend on where the flux, i.e. the integrated current, is measured along the one dimensional 
system, as long as particles cannot accumulate. For large L and J of order j^, Gl satisfies, the 
following scaling 

GL{J;Pa,Pb)~YF{LJ;pa,Pb) (1.2) 
The scaling (jl.2p implies that for large L all the cumulants of Q{t) are of order 1/L, i.e. 

hm ~ -j'^nyP'^^Pb) ■ (1-3) 

Explicit expressions of the Kn{pa,Pb) have been obtained [5l [7] in terms of the diffusion constant 
D{p) and the conductivity (t{p) [23]. One can also show that the large deviation function Gl of 
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the current satistifies the fluctuation theorem [ISlIISlEOlEIlESlIIllIinilSlElE], i.e. 

GL{J;Pa,Pb) - GL{-J;Pa,Pb) = J[^Ogz{pb) - log Z{pa)] 

where z{p) is the fugacity of a reservoir at density p. 



(1.4) 



Reservoir 
Pa 



Reservoir 
Pb 



Figure 1. We are going to consider the distribution of the current flowing through the dashed 
vertical slit of length £ < L. 



In higher dimension, one can study, as in one dimension, the total current flowing through the 
system from one reservoir to the other, but one can also study part of this current. In this paper, we 
consider the SSEP on a square lattice of size L, with periodic boundary conditions in the vertical 
direction and study the current flowing through a vertical slit of length i < L (see figure [I]) . The 
large deviation function Gl,£{J', Pa, Pb), defined as in (II. ip . depends of course on the size i of the 
slit. One reason for considering the fluctuations of this partial current is that in experiments it 
is often only possible to measure the fluctuations of local quantities and not of global quantities 

[iniE]. 

In two dimensions, when £ = L, i.e. when one considers the total current flowing through the 
system, the large deviation function derived from the hydrodynamic theory satisfles for large L and 
J of order 1 a scaling similar to the one dimensional case [4] 

GL,LiJ;Pa,Pb)^F{J;Pa,Pb) (1.5) 

(this would become L'^~'^F{L'^~'^J; pa, Ph) for a cube of size L in dimension d and J of order L'^"^). 
In the present paper we show that G^^^ cannot satisfy the same scaling (jl.Sp as Gl,l and that for 
large L, if one keeps the ratio h = IjL fixed, then for all < < 1 and J of order 1 

GL,Lh{J;Pa,Pb) ^ as L ^ oo . (1.6) 

While, as in (11. 3p . one expects the cumulants of the total flux Q{t) to have a large L limit 

{Q{rr)c , . 

Kn{Pa-,Pb) ■ (1-7) 

in dimension d), we will see by an explicit 



lim 

T^OO 



(which would become ^{Q{t)'"')c — L'^ "^i^niPa^Pb 
calculation of the second cumulant that for (. = Lh, 



lim 

r— >oo 



logL 



as 



L 



oo, 



when < h < 1 , 



(1.^ 



where Q^^^ (r) is the flux of particles through the slit during time r. 
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The fluctuation theorem, which is satisfied as written in (|1.4p for the two-dimensional SSEP 
when J is the total current through the system (i.e. when i = L), has in fact no reason to remain 
valid for i < L: in the large L limit, the difference Gi^ihiJ) — G'l^l/i(— J) vanishes when < h < 1 
so that (II. 4p cannot hold and a singular dependence can be expected in GL,Lh{J) when h ^ 1. 
In Appendix A, we give a simple example of a two site model where one can see clearly that the 
fluctuation theorem is satisfled when one looks at the total current but is no longer valid when one 
considers only part of the current see [l] for a discussion on the validity of the fluctuation theorem 
for partial currents. 

The rest of the paper is organized as follows. In section [21 we recall the hydrodynamic large 
deviation theory [21 [6] and show the asymptotics p.6p . Although the hydrodynamic large deviation 
theory does not predict the correct scaling of the current deviation, the analysis of section [2] suggests 
that local current fluctuations are dominated by vortices. Restarting at the microscopic level, the 
variance of the integrated current is computed for the SSEP on a general graph (section [3|) and 
explicit expressions are obtained for the current through a slit for the SSEP on a two-dimensional 
torus (section H]) . Our exact expression leads to the asymptotics of the form (II. 8p and are compared 
to the results of numerical simulations. Finally the appendices are devoted to comments on the 
fluctuation relation (II. 4p for partial currents, and to some technical calculations. We note that 
sections [21 [3] and Appendix A can be read independently. 

2. Vortices and current fluctuations 

For simplicity, we briefly recall the large deviation hydrodynamic limit theory in the framework 
of the two-dimensional SSEP on the square lattice in the periodic domain A = {1,L}^. At the 
microscopic level, each particle jumps randomly with rate 1 to a nearest neighboring site and the 
jump is allowed only if the neighboring site is empty. After rescaling space by 1/L and time 
by the macroscopic density p{x, t) obeys the diffusion equation |24^ [TS] in the macroscopic 

domain A = [0, 1]^, (with periodic boundary conditions), 

dtp{x,t) = Ap{x,t), x = {xi,X2)eA, (2.1) 

where A denotes the Laplacian. One can also deflne a macroscopic current j (2;, i) = iji{x,t), j2{x,t)) 
in the directions ei , 6*2 which has to satisfy 

dtp{x,t) = -V ■j{x,t). 

The rescaled current j is such that if q(i^i+ea)i'^) is the microscopic integrated current through the 
bond (i, i + e^) (with a = 1 or 2) over the microscopic time interval [0, r], then for a system of size 
L and times r of order L^, one has q[i^i+g^){T) = L J^^^ ja{x,t) dt. 

Using the hydrodynamic large deviation theory, we are going to show that the scaling of the 
large deviations is different for the current flowing through the whole system or through a slit (as 
in flgure[l|). 

2.1. Total current deviations. We denote by Q{t) the integrated total current during the mi- 
croscopic time interval [0, r] through a vertical section of the whole system, say the current flowing 
through the edges {{L/2,i2), {L/2 + l,i2)}i ^L- The corresponding large deviation function 
Gl,l is deflned by 

hm - - log Vto(^^.j]= Gl,l{J) , (2-2) 
T— >oo r \ T J 
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where Pro |^-^^ ~ Jj denotes the probabihty of observing a total current J in the ei direction 

averaged over the microscopic time interval [0,t]. According to the large deviation hydrodynamic 
theory, one expects, in accord with (jl.Sp that lim/,^00 J) = F{J) where the function F{J) = 
liniT^oo Ft{J) with 

W)=mf{ilHi.P)}, and Irij, ,) = \ [ f^Mi. ^ l>..P?Hn^ + l>«p)\ (,,3, 

The minimum is taken over the macroscopic evolutions {p[x,t), j{x,t)} during the macroscopic 
time interval [0, T] which satisfy the constraints 

dtp{x,t) = -V ■ j{x,t), and J =^ dtdx2 h(^(^,x^ . (2.4) 

Remark 2.1. Note that the mathematical statement from the hydrodynamic limit theory [5] relies 
on a more involved asymptotic with a joint space/time scaling: instead of ( f^.^j) . the large deviation 
function for a total current J over the microscopic time interval [0, L^T] is given by 

where Ft has been introduced in 12. 3\) . When writing i2.2\) . i2.3\) . we assumed that in the previous 
expression the limits L ^ oo and T ^ oo can be exchanged. 



As we consider in this section a system with periodic boundary conditions and no sources, the 
steady state is the equilibrium one in which all configurations with a specified total of number 
particles have equal weight. The mean current through the system is therefore and we are going 
to show that for any current deviation J ^ 



F{J) > 0. 

Expanding It in (12. 3p and using Jensen's inequality leads to 



^tU, p) 



1 



T 



dt dx 



Jk 

T 



2p{l - p) 2p(l - p) 
j'^ ^ dtdx {hf + Ct^t(^^ £ j^dt 
where Ct is the contribution of the cross terms in (|2.3p 



+ Ct 



(2.5) 



(2.6) 



dx ji +Ct , 



Ct 



dt dx 



P(l - P) 



dx{S{p{x,Q))-S{p{x,T))} 



with S{p) = — [p \og{p) + (1 — p) log(l — p)]. As it is equivalent to measure the total current through 
any section of the system, the constraint (j2.4p on the current deviations becomes 



J 



J\ 



dt dx ji (x, t) . 



(2.7) 



Thus infj^p |27^(j, p)| ^ TJ^ + Ct- As Ct remains bounded in time, ()2.5p follows from ()2.3p 
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2.2. Partial current deviations. The functional Tt defined in (j2.3|) should in principle provide 
the large deviations of the current through any macroscopic region of the system. We consider 
now a slit of macroscopic height h < 1 (the segment [(1/2, 0), (1/2, /i)]) and denote by Q^^\t) 
the integrated current through the slit during the microscopic time interval [0, r], i.e. the current 
flowing through the edges {{L/2, 12), {L/2 + 1, i2)}i ^ 12 ^ hl- Then, the large deviation function for 
observing a current deviation J 7^ is given by 



lim log Fro 

T— >oo 7" \ T 



J] =GL,LhiJ). 



One expects from (jl.Sp . that limi^oo GL.Lh{J) = Ph{J) with 

F,(J)= hm inf (^lTij,p) 
where ItH^p) is defined in (j2.3p and the macroscopic constraints (j2.4p are replaced by 



(2i 



dtp{x,t) = -V ■ j{x,t), 



and 



We are going to show that in contrast to (|2.5p , the large deviation function Ff^ in 
< /i < 1 (as claimed in p.6p ). 



(2.9) 
vanishes for 



One can bound (j2.8p by 



T 



where the functional Fh{J) is the restriction of Zt to time independent density and current profiles, 



F,(J)=inf i 



dx 



2p(l - p) 2p{l - p) 



where the density and current constraints satisfy 
= V • j = diji{x) + 52j2(x), and 



J 



f 

Jo 



dX2 jl ( 2'^2 



(2.10) 



(2.11) 




Figure 2. On the left, two vortices located at the edges of the dashed slit are depicted. On the 
right, a blow up of one vortex concentrated on the disk of radius (ri,r2): the current (|2.12p . (|2.13p 
is proportional to 1/r for r in (ri,r2) and vanishes outside this annulus. 
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A guess to bound ()2.10p is to consider the equihbrium density (uniformly equal to the constant 
density p) and a current deviation consisting of two vortices localized at the edges of the slit (1/2, 0) 
and (1/2, /i) (see figure [2]) 



VxgA, j{x) = j(^^x-{l/2,h))-^[x-{l/2,{)))), (2.12) 
where $ denotes the vector 

Vx = (XI, X2), a>(x) = ^^^^(^ ^2^^! ( - X2,x0 , (2.13) 

with ri < r2 ^ 1. One can check that the current defined in (j2.12p satisfies the constraint (I2.1ip . 
Furthermore, we can bound Fh{J) by, 

Fh{J) ^ / dx [{nf + {hf] = , ,7 . , , • (2.14) 

4p(l - /)) 7a 4p(l - p) log(r2/ri) 

Letting ri,r2 go to while ^ diverges, we find that Fh[J) = so that the large deviation cost in 
(ES]) is 0. 

Remark 2.2. On a finite lattice of size L, the current has to flow through the bonds and therefore 
the ratio r2/ri is at most L. This imposes a cut-off and the computation (KJ^ based on ^2.10\) 
leads to GL^Lhi J) of order i^^- This logarithmic dependence will be confirmed for the SSEP by a 
direct computation of the current fluctuations in sectionsl^ and\^ 

Expression ()2.14p shows that the cost of the fluctuations due to the vortices is low and one may 
wonder if the vortices are the optimal minimizers of (j2.10p or whether one should expect a more 
complex structure. To check this, we first note that the current j in ()2.1ip is divergence free, thus 
it can be represented as the sum of the curl of a vector field (0, 0, ^(xi, X2)) and a constant vector 
field (Ci,C2) 

j = V X (0, 0, ^) + (Cl, ^2) = (52 -di^) + (Cl, C2) (2.15) 
and the current constraint (j2.1ip becomes 

J = ^(l/2,/i)-'f(l/2,0)+Ci/i. (2.16) 
Choosing the density equal to p, (j2.10p reduces to the varitional principle 

FhiJ) ^ ^ inf^ A L'^^ (^1^)' + (^2^)' + + Cl] . (2.17) 

4p(l - p) *,(Ci,C2) Ua J 

where ^ and Ci satisfy the constraint (j2.16p . The solutions of the above variational problem will 
satisfy 

A^'(x) = a (5x,(l/2,h) - ^x,{l/2,0)) , 

where a is the Lagrange parameter associated to the constraint (|2.16p . This would lead to a ^' 
which diverges logarithmically at the edges of the slit and therefore cannot satisfy condition (j2.16p 
for any non-zero a. Nevertheless, using a cut-off similar to ri in (j2.13p . we can recover the vortex 
like structures (|2.12p . 

Remark 2.3. For more general diffusive systems the hydrodynamic large deviations are governed 
by functionals of the type h2. S^) which depend on diffusion and conductivity matrices [2] . One could 
extend the previous discussion to these cases and the large deviation function Fh{J) would vanish 
as soon as h < 1. For open systems, similar computations can be done as the fluctuations are 
dominated by vortices localized at the edges of the slit /12.12\} and the reservoirs play a negligible 
role. 
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Remark 2.4. We note that in analogy to ^2.7\) we can consider the partial current specified in 



as a limiting case of an integrated current in a domain B C A, 



1 



T 



JB = 7f / ji{x,t)dx. (2.18) 
^ Jo Jb 

Taking B to be a rectangle of height h and width w we get that the flux through the line segment h 
is given by w~^Jb, in the limit — > 0. For the large deviation of Jb one can repeat the analysis 
leading to i2. 6]) yielding 

J2 

inf Jt(j, p:B)^ T-§- + Ct ■ (2.19) 

3,P \B\ 

This non-zero lower bound reflects the fact that any vortex flow used to implement the fluctuation 
Jb will have to be of a size w or greater. 

3. Current fluctuations on a general Graph 

In this section, we consider the SSEP on a general connected graph {A,£/^) where A is a finite 
set of sites and £\ the set of edges on which particles jump with rate 1 according to the exclusion 
rule. We also suppose that particles are created and annihilated at the sites in the subset of A 
(Q may be empty). For any site i in 0, we suppose that creation and annihilation occur at rate Oi 
(3i (and for simplicity we choose + /3j = 1). In section [H we will apply the results obtained for 
general graphs to the microscopic domain A = {1, L}^ and derive explicit expressions in this case. 

If {i,j) is an edge in then the number of particle jumps from i to j during the time interval r 
is denoted by (?(jj)(T). The current flowing through (i, j) during time r is then g(jj)(T) — (?(j^j)(T). 
If creation and annihilation occur, we enlarge the graph {A,£\) by associating to each site i in 
a new site i. The site i can be interpreted as a source and we denote by q(j^i-){T) the number of 
particles created at i and q(^i~i){T) the number of particles annihilated at site i. It is convenient to 
enlarge the graph into i?a by adding to the original graph the new edges {i,i)i^n and {i,i)i£n. 
We denote by the set of the new sites, and by A the union of A and 0. 

For any field {^fe}?) indexed by the edges and such that for any edge {i,j), = — w4(j j), we 

are going to study the fluctuation of the integrated current defined by 

Q-^(r)= 5] A%(r), (3.1) 

where the sum is over all the oriented edges b. The field Ab can be thought as a test function. 

One can define the divergence and the gradient on the graph. For any field and any site 

i in A 



dwAii) = Y,Ai,j)^ (3-2) 



where the sum is over all the edges leaving site i (this includes the edges {i,i) if creation or 
annihilation occur at site i). For any function Hi in A, the gradient is a function indexed by the 
edges b = {i,j) in £a 

VbH = Hj - Hi . (3.3) 
In the following, we will consider only functions if in A equal to in 
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3.1. Gauge invariance. Before, computing the variance of Q'^(r) defined in (f3lD . we first show 
that for large r similar asymptotics of log(exp (AQ-^(r))) can be obtained for different choices of 
{Ab}b- 



For any site i in A, one has with notation (13. 2p 

rj{i, t) - rj{i, 0) = ^ - W • (3-4) 



jeA, 



This implies that for any function if on A (equal to on Q) 

\r)=^VbHq,iT)=^H, f ^ ^O.^W - (?(i,,)(r) I = J] i/, (r?(i, r) - 0)) . (3.5) 

jeA, / ieA 



fcefA «6A 



Therefore, Q^^{t) remains bounded when the time diverges. As a consequence, for any A G M, 
and Hj 



lirn^ - log / exp I A ^ ^^(r) j \ = lirn^ i log / exp | A ^ (A + V^ii) qbir) j 



(3.6) 



Thus the large deviations of Q"^(t) and Q"^"'"^^(r) with respect to time are the same. 

We are going now to recall how a field Ah in £\ can be decomposed as 

Ab = Vb + VbH, (3.7) 

where is a function in A (equal to on Vl) and V is divergence free in A 

\/i £ A, divV{i) = 0. 

Note that no conditions are imposed on divy(i) for i in J7, if there are sources (il ^ 0). If there 
are no sources (fi = 0), then H is defined up to a constant. 
For decomposition l\3.7\i to hold, H has to be the solution of 

Vi G A, div^(i) = AHi = ^{Hj - Hi) . (3.8) 

In the case with sources {il. ^ 0), the solution of ()3.8p can be written in terms of the Green's 
functions, defined for any site k by 

ViGA, AGf^ = -(5i,fc, and Vi G Gf^=0. (3.9) 

Thus for j G A 

Hj = -Y^ divAik) Cf^ , (3.10) 
fceA 

and the field Vb = Ab — V^H is divergence free. 
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3.2. Variance of the current. We are going to compute the variance of Q'^(r) = X^;,g£a ^Ibi''')- 
Prom (j3.6H3.7p , we know that to compute large time asymptotics it is enough to consider A which 
is divergence free. 

One has 

a,log(exp (AQ-4(t))> = 

X- / nA ^ .\ ((l-r?,)exp(AQ-^(r))) ^ ^ ^ (r?, exp (AQ-^(r)) ) 

|;a.(exp(A^(.,)) - 1) ^exp(AQ^(r))) + (3^{eM-XA^l,)) " l) (exp (AQ^(r))) 

^ V- . ^ ,, (r?i(l-r?,)exp (AQ-^(r))) 

+ 5: (exp(A^(,,)) - 1) ^,,^^^^A^,^^^ ' (3-11) 

(«j)G<?A 

where the sum is over all the oriented bonds and (•) denotes the average over the random 
process in the time interval [0, r] and over an initial condition chosen according to the invariant 
measure for the SSEP. The procedure to derive (|3.1ip is similar to what was done in [11]. One 
considers all the possible moves occurring during an infinitesimal time interval dr and their contri- 
butions to (exp (AQ"^(r))^. The first terms in (|3.1ip correspond to a jump of a particle from site 
i to site i (creation) or from site i to site i (annihilation), whereas the last term corresponds to a 
jump from site i to j. 

Let us denote by {■)x the expectation of the tilted measure: for any function / 

^^^^^^^ - ™ (exp(AQ-^(r))) " ^^'^^^ 

Using the symmetry ^(jj) = — w4(j j) and the relation + /Jj = 1, we get by expanding ()3.1ip 
for small A 

hm dr log (exp (AQ-^(r)) > = A J] A^j^^^ai + A ^ div^(i) (r/^A (3.13) 
For A small, one expects that 

{v^)x = {m) + o(A), - n,))x = - m)) + o(A), (3.i4) 

In principle one would need to know the first order correction to {rfi)-^ to obtain (13. lip at the second 
order in A. 

For A divergence free, the term {rji)\ disappears and the formula (j3.13p simplifies 

1 A^ 
lim - log (exp (AQ-^(t))) = AMean + — Var + 0{\^) . (3.15) 

T— >oo T 2 



where 



r— ►oo 



Mean = lim S^±±Jl = An i 



Var = ^hm S^^JUJl = ^ - Vj)) + E(Alo)' (^^(1 " + <^ " ^^)) '(3-16) 
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where the sum is over all the oriented edges. If there are no sources then the second term in ()3.16p 
disappears. 

For a general field A, we can use the decomposition ()3.7p . If 7^ 0, there is a representation of 
V in terms of Green functions (j3.10p : for any S A 



(3.17) 



fceA 



Invariance (j3.6p implies that the asymptotic formula for (|3.17p is given by ()3.15p with A replaced 
hyV 



Mean = ^ V^^ .^a^ 



iefi 

(j,j)e£-A ien 
where as before the sum is over all the oriented edges. 

Further simplications can be obtained if the system is in equilibrium at density p, i.e. if the 
intensities of the sources are such that Oi = p, (3i = 1 — p. In this case, (?7i(l — r/j)) = p(l — p) and 
by expanding p.lSp with Vb = Ab — ^bH, we get 



Var = p{l — p) 



,(*j)GfA («j)e£'A (*j)GfA 



E 4,,)+4Ediv^»^^-2E^^*^i 

.(«j)G^A »eA ieA 



where the second equation is obtained by summation by parts. From (j3.8p . one has Aif,, = div^(i) 
so that 



Var = ^(l-p) 4,,)+2j]div^(z)i7, 
Replacing H by ()3.10p . we finally obtain 



(3.19) 



Var = p{l — p) 



E 4j)-2E^J'^div^(^)div^(i) , (3.20) 

_(ij)G£A i.fceA 

where the sum is over all the oriented edges (i,j). 

4. Two DIMENSIONAL SSEP 

In this section we will apply the general results of section 3 to the SSEP in the periodic square 
lattice A = {1,L}^ with nearest neighbor jumps and derive explicit expressions in this case. We 
consider the integrated current flowing through the edges in F^ = {L/2, L/2 + 1} x {1, ^} given by 



{j,i+ei)er^ 



(4.1) 
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with 6*1 = (1,0). The integrated current can be rewritten as Q"^ defined in with 

fi, if (i,j) = (i,i + ei) eri, 

Vi,j e A, A,i = < -1, if («, j) = (^,^ - ei) e 
I 0, otherwise . 



L ' 



(4.2) 



The gauge invariance (13. 6p is easily illustrated in the two dimensional case. Let and be 
the 2 sites of the dual lattice such that F^^ is the set of edges intersected by the segment (z^, z~[) 
(see figure [3]). Let 7 be another path connecting z'l to z~^ on the dual lattice, then we can define 
the current Q'^(r) flowing through the edges intersecting 7, where B generalizes ()4.2p for the edges 
intersecting 7. One can check that 

A = B + VH, (4.3) 

for some H. Therefore, the statistics of the currents Q"^(r) and Q'^(t) are asymptotically the same 
at large times (|3.6p . 

z+ 7 
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Figure 3. The dashed hnes represent the dual lattice and is depicted by the thick edges. The 
function H defined in equation (|4.3p is equal to 1 in the grey region and outside. 



4.1. Computation of the variance. We turn now to the computation of 



Var 



lim ■ 

r— »oo 



the asymptotic of the variance of Q^{t) = Q"^(r) (j4.2p for large r. On the periodic domain, 
the variance is given by (j3.18p without the source term. As the invariant measure is uniformly 
distributed, (r?i(l — rjj)) depends only on the number N of particles and the size L. Let Sl^n = 

'l^Il^-V) — {Vi{^ ~ Vj)) foi' ^ / then the expression (j3.19p remains valid 



where H is given by (13. 8p which reads now 



(4.4) 



1 < s < . 



Vi G A, AHi = div^(i) 



ifiGr/ = {L/2,s} 
if i G r^- = {L/2 + ei,s} 
otherwise 



1 < s < . 



(4.5) 
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Thus H is equal to 



Hi 



ik,k+ei) 



fcGF 



(4.6) 



where the Green's function (j3.9p is replaced for any sites k, k' in A by 



Prom (|4.4|) . we finally obtain 



(4.7) 



Var = 2Sl. 



N 



E 



^{k,k+ei) ^{k,k+ei) 



(4.8) 



The Green's function (|4.7p is given for any i = (ii,i2) in A by 

^(fc,A:') _ J_ ^ exp (i g • (i - A;)) - exp (i g • (z - /c')) 



L2 



4 — 2cos(gi) — 2cos(g2) 



(4.9) 



gi,<727^(0,0) 

where q ■ j = Qiji + (72^2 stands for the scalar product with qi = 27r^,q2 = ^vr^ for mi,m2 in 

-cos (92) 



{0, L — 1}. Thus (j4.8|) becomes with the convention that \ = fo'^ 92 = 



Var = 2Sl. 



N 



(1 - COs(gi))(l - COs(g2^)) 



91,9^^0,0) " ^°s(92)) (2 - cos((7i) - cos((/2)) 



(4.10) 



where the subscript has been added to keep track of the dependence in L, L Using the identity 

l_ ^ 1 - cos(g2^) _}_y^ 1 - cos (2^£^) ^ _ ^ 
L ^ 1 - cosfoo) ~ I 1 ~ ~ ^ 



92 



m=l 



COS 



we finally rewrite (j4.10p as 



e 1 



1 - cos {q2l) 



Var = 25,,^!^ + ^ E 2-cos(,,)-cos(,2) 

91,927^(0,0) 



(4.11) 



with qi = 2TT^,q2 = 27r^ for mi, m2 in {0, L - 1}. 



One can show, see Appendix B, that for large L and i, expression (j4.1ip becomes for h = i/L 



Var 



2Sl. 



N 



IT 



log L + + loe 



sinh(7r/i) \ 3 log 2 



TT 



+ 



+ 7 + ^ log 1 + 



sm 



(nh) 



m > 1 



sinh (vrm) 



, (4.12) 



where 7 ~ 0.577 is the Euler constant. 
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4.2. Time dependence. To compare with the results of simulations, it is useful to calculate how 
the moments of Q'^(t) depend on r. At finite times r, the gauge invariance ()3.6p is of no use. We 
focus now on systems with no sources and study the variance of Q"^(r) for a general field A at 
finite time r. 

Following the same procedure which led to (IS.llh . (I3.13p . we get up to the second order in A, 
dr <exp (AQ-^(r))> = A ^^dwAii) (r/i(r)exp (AQ-^(r))> 

+Y E (Am-))'(^^W(i-^.-W)«^p(^'Q-^W)> ' (4-13) 

where the second sum is over all the oriented edges. We therefore need to determine {rn (r) exp (AQ'^(r)) ) 
to first order in A. To do so, we can write 

dr (r/.(T)exp (AQ-^(t))> = (exp(A^(fc,,)) - 1) (r?,(T)7?fc(T)(l - r?,(r))exp (AQ^(r))> 

('=j)s£a 

+ ^exp(A^(,,)) (r/,(r)(l - r?,(r))exp {XQ^{t))) - (r/,(r)(l - r?,(r))exp {XQ^{t))) , 

where the first sum is over all the oriented edges {k,j) which do not intersect i, but the second sum 
is over the neighbors j of i. As for (jS.lip . this expression can be derived by adding the contributions 
of all the single moves which may occur during the infinitesimal time interval dr. Expanding to 
first order in A, we get for a given r 

dr{vi{T)exp{Xq-^{T))) = X Akd){^ii^k-Vj)) (4.14) 

('=j)e£A 

+ Yl hir)exp (AQ^(r))> - (r7,(r)exp (AQ^(r))) + A^(,.,) (7?,(1 - ry,)) . 

In the periodic case, the first term in (j4.14p vanishes since — rji)) = Sl,n, is independent of 
Hence 

dr (r?.(r)exp (AQ-4(r))) = ^ (r?,(r)exp (AQ-^(r))> - (r/i(r)exp (AQ^(r))> 

i; (jj)G£A 

+XSL,Nj2Aj,i)- (4-15) 

We introduce for any site A; in A the time dependent Green's function, solution of 

Vi G A, drG^^} = AgI'} + = (g15 - ai'J) + 6,,k , (4.16) 

with the initial condition g'^qJ = 0. Integrating (j4.15p . one obtains to first order in A 

(7?,(r) exp (AQ-^(r)) > = + XSl,n ^ E Am = im) - ^Sl,n ^ Gi'}divA{k) . (4.17) 

k j~k k 
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Using (liTH) in (fiJ3]) . we get 

5,(exp (AQ^(t))> = A^Sl, 



for the second order term in A 

\ E (AM))'-E<idiv^(^)div^(.)] , (4.18) 



.N 



where the first sum is over all the oriented edges. 

For the periodic square lattice A = {1,L}^, the Green's function ()4.16p is given for any k = 
{ki,k2) and i = {11,12) in A 

1 1 _ -(4-2cos(gi)-2cos(g2))T 

^ [ ,u<^^m 4-2cos(gi)-2cosfe) 

with qi = 2Tr^, q2 = for mi,m2 in {0,L — 1}. Using (j4.18p . we deduce an exact expression 

for the variance of the current Q^{t) through a slit {A is given by (j4.2p ) 

mrf). = 2S,.J\sL-' E (l-e-''---''-'--"';'»-)(l-cos(..))(l-coste.)) ' 
Jq \ 2 - cos(gi) - cos(g2) l-cosfe) 



We finally get 

Mll!)£ = Var + ^ y (1 - e-(4-^^°^('?^)-^^°^fa^))") (1 - cos(gi))(l - cosfe^)) ^^^^^ 
T tL"^ ^ (2 - cos(gi) - cosfe))^ 1 - cosfe) 

,92^(0,0) 

where Var, given in (j4.10p or (j4.1ip . is related to the large r asymptotics. 




Figure 4. ' is measured versus I in numerical simulations for the SSEP on a square of 80 x 80 

sites for times r = 250, 750, 2500 (the results decrease with r). The continuous lines represent the 
theoretical predictions (|4.19|) at these times, as well as the limit r = 00 (|4.11|) . Expression (|4.19p 
fully agrees with the simulations. 



4.3. Numerical simulations. We show now the results of the simulations of the SSEP on a square 
lattice of size L = 80 with periodic boundary conditions at density /> = 1/4 (without reservoirs). 
The initial condition is chosen at equilibrium (i.e. the Lp'p particles are put at random positions on 
the square lattice). For each simulation, we measured the flux Q^{t) through a slit of microscopic 
length i during time r (see figure [Hand ()4.ip ). 
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In figure [U our data for ^ ' are compared with the predictions obtained from (|4.19p for 
different times r = 250, 750, 2500. The simulations are averaged over 10^ realizations. We see that 
unless the time is long enough, the results differ significantly from their infinite time limit (j4.1ip . 
One can notice that for short times, the variance grows essentially linearly wrt £ as the current 
fluctuations are simply the sum of the (almost) independent contributions of the local current 
fluctuations along the slit. 

In figure the theoretical curve Var = limT-^oo ^'^'^ ^ computed in (|4.11|) is shown for several 
system sizes L = 40, 80, 160, 320. One can notes that the variance of the current flowing through 
the whole system {£ = L) is independent of L. 

1-2 I 1 1 1 1 1 1 1 



1 - 




0.2 - 

I I I I I I I 

50 100 150 200 250 300 350 



Figure 5. Theoretical prediction (|4TT1l of lim^_oo ^ ' ' versus ^ for L = 40, 80, 160, 320. 

In figure El the same data as in figure [5] are shown but the horizontal axis is now l/L. One can 
see that for large L, Var grows linearly with logL as predicted in (j4.12p . 



1 - 




0.2 0.4 0.6 0.8 1 

l/L 

Figure 6. Theoretical prediction of linir^oo ^'^^^^'''^ versus t/L for L = 40, 80, 160, 320. 

5. Conclusion 

In this paper, we have computed the variance of the local current for the symmetric simple 
exclusion process on general graphs with reservoirs (13.180 or without ()3.20p . In two dimensions, 
our exact expression leads to the asymptotics of the variance through a slit ()4.12p . The logarithmic 
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dependence of the variance confirms that vortices dominate the local current fluctuations. As a 
consequence a fluctuation of the partial current, say the current flowing through 99% of the system, 
does not obey the same scaling as a fluctuation of the total current. For two dimensional diffusive 
models, we have also seen that the hydrodynamic large deviation theory does not catch the correct 
scaling of the current deviations ()1.6p . Finite time corrections to the variance were also computed 
(I4.19P and compared to numerical data (figure H]). Finally the fluctuation relation ()1.4I) for partial 
currents is discussed in Appendix A. 

It would be interesting to investigate the scaling of partial current deviations in higher dimen- 
sions. Another challenging issue is the computation of the full large deviation functional for partial 
currents. 

Appendix A: The fluctuation theorem and partial currents 
For the total current, the fluctuation theorem (jl.4p holds (in any dimension) 

GL,LiJ;Pa,Pb) - GL,L{-J;Pa,Pb) = J[log z{pb) - log z{pa)] ■ (A.l) 

This fluctuation relation is based on a global symmetry: the fluctuation to produce the current — J 
is simply the time reversal of the fluctuation to produce the current J. One may wonder how this 
generalizes to the function G^^i. In this appendix, we show by considering a very simplified model 
that the fluctuation relation (jA.ip is in general not satisfied for partial current deviations. 

We consider the SSEP with two sites {1,2} connected to reservoirs. At site 1, creation (resp 
annihilation) occurs at rate a (resp 7) and at site 2, creation (resp annihilation) occurs at rate 5 
(resp (3) (see figure [7]). The exchanges between sites {1,2} obey the usual exclusion rule, but they 
can occur through two edges with rate 1. On the one hand the model behaves like a SSEP with 
exchange rate 2. Thus the total current flowing from site 1 to site 2 obeys the fluctuation relation 
(lA.ip (with Pa = and pb = sTfi^' other hand one can also consider a current deviation 

J through one of the two edges. Heuristically one can see that the system is going to use different 
strategies to produce a current J or — J. Imagine the extreme case with only creation at site 1 
and annihilation at site 2 (7 = and 5 = 0). For the total current, there is no way of producing 
a negative flux and the relation (jA.ip is degenerate: log z{pi,) — log z{pa) = —00. On the other 
hand, a negative current can be achieved through the lower edge by letting a single particle cross 
the lower edge from site 2 to site 1 and then use the upper edge to go back to site 2. This latter 
mechanism mimics the vortices discussed in section [2l Thus, the fluctuations to produce current 
deviations J or —J are not related by time reversal. In general, both mechanisms (total current 
deviation and local vortices) combine and there is no reason to expect a symmetry such as (jA.ip . 




Figure 7. A reservoir at density pa (resp pb) is acting on the left (resp right) site by creating 
particles at rate a (resp S) and annihilating particles at rate 7 (resp j3). We consider the large 
deviations of the current Q'^ flowing through the lower edge. 
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We analyze now the toy model analytically. We define Q'^ as the integrated current flowing 
through the lower edge during the time interval [0,t] (see figure [7]). Instead of trying to check an 
expression like (|A.ip for the large deviation function we look for a symmetry at the level of its 
Legendre transform. As in [TT], one knows that 



\ 



VA, lim - log (exp {\Q'^ > = /x(A) , 

where /i(A) is the largest eigenvalue of the operator 

( -a-5 13 

a -7-5-2 l + e~^ /3 

<5 1 + -a-/3-27 

\0 5 a -1- fi ) 

The fluctuation relation (jA.ip would say that there exists a constant E such that 

VA, ^(-A - 2E) = i2{X) (A.2) 

In order to prove that the previous relation does not hold, we consider for simplicity the case 
a = 2, 7 = 1 and 5 = 1, (3 = 2. Then the characteristic polynomial of Lx is 

P{u) = {3 + u) {16 + u{U + u{13 + u)) - 2(8 + u)cosh[A] - 6sinh[A]) 

For ()A.2p to be satisfied, /u(A) should be a root of this polynomial and of the polynomial associated 
to This implies that 

2(8 + /x(A)) cosh[A] + 6 sinh[A] = 2(8 + /i(A)) cosh[-A - 2E] + 6 sinh[-A - 2E] 

leading to 

VA, ^(A) = -3coth[S] - 8 or /i(A) = -3. 

As /i(A) cannot be independent of A, we obtained a contradiction. Thus the fluctuation relation 
()A.2p does not hold in this toy model. 

For the total current a similar calculation shows that /^(A) is the root of 

Q{u) = (3 + u){20 + u{6 + u){7 + u) - 20cosh[A] - 12sinh[A]). 

This is invariant under the symmetry A — > —2 log 2 — A, implying that (|A.2p is satisfied. 

Appendix B: 

In this appendix we derive the large L,i expression ()4.12p of the variance (|4.1ip . 
Define In and Jn by 



TV 



1 



n=l 



ri' 



+ 62 



and 



Jn 



N 



dx 



62 + x2 



One has for large N 
In = Jn + 71"- 
Recall also that 



1 



1 1 

27r6tanh(7r6) ~ 2^ ~ 2^ 

dx 



+ o(l) = J;v + TT^ 



7r6[exp(27r6) - 1] 27r262 



2-B- cos(27rx) ^{1-B){3-B) ' 



+ o(lX5.1) 



(5.2) 
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From ()5.1U5.2p . one can show, by taking 6^ = L^(l - cos q'2)/(2vr^), that for < ^2 < 27r 

1 1 



1 



^ nf^O 2 ~ ^"^^2 - cos ^ v^(l-COS g2)(3-COS ^2) 



1 



sin(f )(exp[2Lsin(f )] - 1) 



+ o(l). (5.3) 



The main contribution to the difference between the sum in ()5.3p and integral ()5.2p is given by the 
terms with ni close to or to L. In both cases (1 — cos ) can approximated by its second order 
expansion and the last term in (j5.3p is obtained thanks to (|5.ip . 

can be rewritten as 

1 1- 



3— cos 92 



- + - 



i V - 

^„fr'o2-cosg2 -cos^ ^(l-cosg2)(3-cosg2) 2sinf sin(f )(exp[2L sin(f )] - 1 
One can then perform the sum over q2- For large L, the first term becomes an integral 



-+o(l) 



/ dx- 

Jo 



1 



3~cos(27rx) 
2 



1 



1 - \/l + sin^ , 



10 y/[l - cos(27r2;)][3 - cos(27r3;)] 27r 

For large L one can also show that 

L-l 



sm 



(/> \/l + sin^ 



log 2 
27r 



Z ^ 2sin^^ ^ ^ 



n=l L 

For large Z and L with Z = L/i 



log L + log ( - ) + 



n=l " L 

There is also the identity 



1 cosf— Z) 1 1 
7 E 9 - ""nn ^ log (2 - 2 cos(27r/i)) = — log(2 siniTTh)) 
L ^ — ^ z sm -f^ zvr vr 



E 



2[1 - cos(27r7i/i)] 1 



^ n7r(e2"'^ - 1) 

Putting everything together one gets that 
^ L-l L-l 
Yl 



m=l 



sin^(7r/i) 
sinli^fvrm) 



1 - cos{q2l) 



L2 

1 

TT 



2 - cos 0^2 - COS qi 

ni=0 712=1 



log L + loe 



sin(7r/i)\ 3 log 2 



TT 



+ 



+ 7i? + log ( 1 + 



sin^ (vr/i) 



, sinh (vrm) 

m ^ 1 ^ ' ' 

Note that for h small, i.e for 1 <^ I <^ L, one recovers a well known expression (see [17] page 198). 
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